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General Introduction 

This is the Tirst annual technical report -for research carried 
out under NASA Contract NASW-3477 by J. S. Lewis Associates, Inc., 
a not— For-prof i t corporation chartered in the Commonwealth oi 
Massachusetts. The Corporation was founded for the purpose of 
carrying out basic research under contract to governmental and 
industrial clients, and is structured to carry out theoretical and 
computer modeling research using minicomputer technology rather 
than large mainframe computers- The Corporation will make every 
effort to keep overhead expenses as low as possible, and will 
encourage recent PhDs with excellent qualifications who are unable 
to find tenure-track university positions in the present depressed 
academic job market to seek research support under its aegis- 

Our work during the first year of this contract has concerned 
three research areas of fundamental interest to the NASA planetary 
exploration program- These areas are the chemistry and physics • 
of formation of the satellite system of Jupiter, the comparative 
study of the conditions of origin of Venus and Earth, and the 
condensation behavior of presolar solids in physically realistic 
as trophysi Cc 1 settings- Our progress in these areas is described 
in sections I, II and III of the present report, respectively. 
The accomp 1 i shmen hs of the first year have closely approximated 
our oriqinal intent as expressed in last yearns proposal: however, 
the means by which these ends were achieved have changed- It was 
our oriqinal proposal to rent a mi ni compute at the outset of this 
work, choosing the most capable one for the available funds. ♦ 
However , ah the time of the outset of this work it became known 
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that a new IBM system (the Personal Computer) would be on the 
market in a few months. Accordingly, we chose to begin this 
project by renting a terminal and buying time on a large computer. 
During the year the IBM computer came out, and we are now running 
on it in Basic, pending the availability and assessment ai a 
FORTRAN compiler. We are delighted with the speed and capacity of 
this computer, and have nearly finished the task of program 
conversion, ’ 

Our cost experience during the first year reflects the nature 
of a newly begun program; secretarial expenditures have run low 
and late in the year because the burden of report preparation 
could only be begun once sufficient research had been done to 
write up. Since we have to date only submitted one report for 
publication, our -publication expenses have so far been virtually 
nil. The budgeted travel money could not be used due to the 
untimely illness of the P. I. None of these circumstances should 

apply during the coming year, and the budget has been planned 
accordingly. 
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SECTION I 


CHEMICAL AND PHYSICAL MODELS OF THE 
JOVIAN SUBNEBULA 


John S. Lewis 

J. S. Lewis Associates, Inc. 
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Abstract 

A semi emp X r i cal physical moct^l o-f the Jovian SLibnebula has been 
developed by analogy with the primitive solar nebula itself. This 
model employs simple po/jer^-law scaling of the dependence of 
temperature, pressure, eensity and surface density with 
jovicentric distance. The chemical aspects of this model are 
developed according to the principles developed in the study of 
the thermochemi stry and gas kinetic behavior of the solar nebula,* 
but with important modifications to take into account the higher 
pressures and densities in the jovian subnebula. These 
differences in . physical conditions are reflected both in th^ 
higher condensation temperatures of pure substances, and in the 
degree to which gas phase and heterogeneously catalysed g^s 
reactions approach equilibrium. In the jovian subnehula, unlike 
the solar nebula, the reduction of the hi gh-temper ature gases CO 
and N2 by reaction with molecular hydrogen can readily precede 
toward equilibrium, thus making both NH3 and CH4 relatively more 
important components in the jovian system. The bulk compositions 
and densities of the inner satellites of Jupiter are calculated 
from this physico-chemical model and compared to observations. 
The recent petrological model of lo which attributes its active 
sulfur maqmatism and sul fur dioxide venting to an initial state of 
rather high degree of oxidation is found to be reasonably 
concordant with this model. It is proposed that Europa differs 
from lo chiefly in that it has suffered a less severe thermal, 
history. The general features of this model are applicable with 
minor modification to the systems of Saturn and Uranus. 
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Introduction 

Advancing theoretical understanding of the early history of the 
jovxan system and in our observational data base for Jupiter and 
its satellites permits an attempt to use the geochemical modeling 
techniques which have been developed to treat the solar nebula in 
this narrower context. Perhaps the greatest shortcoming of any 
model for the early solar system is that it almost always 
contains certain ad hoc postulates or adjustable parameters which, 
however plausible and successful, may be utterly lacking in 
generality. Because of the nature of the observational evidence, 
we may not readily turn to another solar system to test our‘ 
theory. Thus we are left quite uncertain whether we have a 
general scientific theory, or a mere engineering model of our own. 
system which reflects many features we observe, but may lack 

essential physical insights into the fundamental processes at 
work. 


The existence of regular satellite systems about Jupiter, Saturn 

and Uranus provides us with three new opportunities to expose our 
theories to tes:. 


We shall here explore 
between the formation of 
nebula and the formation 
primitive jovi an subnebul a 
applicable to the sysi terns 


a very particular and simple analogy 
the planets out of a primitive solar 
of Jupiter ^5 satellites out of a 
• While this model is in . principle • 
of Saturn and Uranus -- . 


as wel 1 , 


our . 
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greatly superior knowledge about the jovian system compels us to 
test the model -first in that context. 

Following Cameron <1978) and others, we shall assume solar 
elemental abundances for the jovian subnebula. By analogy to the 
mass distribution in the solar nebula, computed by distributing 
the mass of each planet over a torus centered about the orbit of 
that planet, and replenishing missing volatiles up to solar 
relative abundances (Wei denschi 1 1 i ng, 1977), we find a surface 
density, which is proportional to the heliocentric distance, 

r, raised to the -1.6 power. From the general dependence of 
planetary density and bulk composition on hel i ocentr ic distance, 
we find that the temperature in the solar .-i.'bula varied as r to 
the -1.1 power (Lewis, 1974). These assumptions combine to 
require that the midplane nebular pressure, Pc, drops off as the 
-3.6 power of the heliocentric distance, while the midplane gas 
density, drops as the -2.5 power of r. 

The densities of the Galilean satellites drop off monotoni cal 1 y 
with distance from Jupiter, in a fashion reminiscent of the 
behavior of the densities of the planets as a function of their 
distance from the sun- The masses and densities of the Galilec.n 
1 1 tes are given in fable i- By eamparison with the 
theor et i cal 1 V computed bulk density of the condensed portion of a 
solar-composi tioh - system (Lewis, 1972, 1974), it . is plainly- 

evident that water ice must be a major component of both Ganymede 
and Caliisto, but not lo and Furopa. Accordingly, we shai 1 place 
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Table I 


Radi i , 

Masses and Densities 

o-f the Galilean 

Satel 1 i tes 

Satel 1 i te 

Radius 

24 

Mass (10 gm) 

Density (< 

lo 

1Q20 

89. 1 

“T CT'?' 




w . wJO 

Europa 

1552 

48.7 

3.03 

Ganymede 

2635 

149.0 

1.93 

Cal 1 i sto 

2420 

106.5 

1.83 
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the condensation threshhold of ice halfway between the orbits of 
Europa and Ganymede. Now, using the nebular dependence of 
surface density and temperature on distance found above, and 
requiring that the surface density must be sufficiently high to 
provide the observed mass of the most massive Galilean satellite, 
Ganymede, we may then calculate the expected compositions, masses 
and densities of the other Galilean satellites and Amalthea (J5) . 

We assume for the sake of simplicity that condensed matter is 
^•^creted into that satellite whose orbit is closest to the point 
of the solids; this assumption should be very good for^ 
the purpose of estimating the total accreted mass, but less 
adequate for describing the compositional diversity of each 
satellite (and hence its exact bulk density and volatile 
content) . 

Our procedure shall be to calculate first the temperature, 
pressure, and density profiles in the subnebula, and then to 
superimpose on this structure the condensation points of the 
rock-forming and ice-forming species in the nebula. We then 
calculate the bulk composition of each satellite, its mass, and 
its density. We then shall compare these satellite models to • 
observational evidence. Finally, we shall critiq-.e the model 
both in light of the available evidence and the theoretical 
naivete of its as'sumpt i ons. 
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Physical and Chemical Structure the Subnebula 

P-->r results for the temperature profile in the jovian subnebula 
are given in Fig. 1 . The temperatures at the present orbits of 
the inner jovian satellites are found to be about HOOK for 
Amalthea (at a pressure of about 80 bars), 450K for lo (4b), 275K 
for Europa (0.8b), 170K for Ganymede (0.16b), and 140K for 
Callisto (.05b). The ice condensation point along this adiabat 
lies at 225K, ammonium hydrosulfide condenses at 170K, and ammonia 
monohydrate at 150K. The temperature given for Callisto is of 
course uncertain to some deg-ee due to the uncertainty in 
determining exactly where the adiabatic thermal structure of the- 
Inner part of the subnebula gives way to the radiatively 
controlled and nearly isothermal outer portion of the subnebula. 

Thus it is far from certain tnat ammonia monohydr ate would in 
fact condense and be accreted into Callisto. 

While these results were being prepared for publication, a 
preprint on the same subject (Lunine and Stevenson, 1901) was sent 
out. Although the physical model of the nebula emploved by Lunine 
and Stevenson was constructed to quite different criteria than the 
present model, the results are strikingly similar. In general, 
the temperature and pressure gradients in the present model are 
slightly steeper, while the temperatures and pres'>ures of 
formation of each satellite are slightly low-r. 


uince Lunine and Stevenson ha 


ve concentrated on the composition 


relationships between Ganymede and Callisto, and " si 


nee our 
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Caption to Fig. 1 

Temperatures and Condensation Fronts in the Jovian Subnebula 

The vertical solid line marks the present radius of Jupiter, 
and the vertical dashed lines mark the orbital semimajor axes of 
the five innermost large satellites, J5 (Amalthea), J1 do), J2 
<Europa), J3 (Ganyn.ede) , and J4 (Callisto). The condensation 
points for refractory oxides, liquid metallic Fe-l-’i , enstatite. 
albite, troilite, the endpoint for oxidation of fe, and the • . 
condensation points of water ice, ammonium nydrosul f ide, and the 
solid monohydrate of ammonia are indicated by the horizontal 
dashed lines. Condensation behavior is predicated upon the 
assumption of adiabatic structure and solar elemental abundances. 
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results for these bodies, calculated for a nebular mass about 
half that used by them and for a sliqhtly lower temperature, are 
striJang,. similar, we shall concentrate on the predictions of 
this jnodfl -for Aoialthea, lo and Europa. 


Our calculated condensation point for Amalthea lies well below 
the condensation curves of metal, magnesium silicates, and alkali 
aluminosilicates, but well above the formation temperature of 
troilite, FeS.. A very highly reduced assemblage would be 
expected, in which kamacite and taenite would be abundant, but the 
FeO content of the principal silicate lenstatite) would be 
neglKgible. The relationships between the formation condi ti ons ’.of 
the inner satellites and the condensation curves of the 
rock-forming elements can be seen in Fig. 2. Note that the 
composition of Amalthea is predicted to be closely similar to that 
of V-nus, even though the pressure at the point of formation of 
Amalthea is approximately times higher. 


The composition of lo is expected to be closely similar to that 
of Mars, although with a higher initial water content due to the 
formation of the material of lo within the stability field of 
serpentine. Whether serpentine can be formed in the solar nebula 
by reactions between water vapor and f erromagnesi an minerals is 
certainly debatable: it is more plausible to attribute the 
meteoritic occurrences of this mineral to secondary alteration 
reactions which, took place on the meteorite parent body. The 
water vapor partial pressure in the jovian subnebula at the point 
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Caption to Fig. 2 

Pressure-Temperature Structure of the Jovian and Solar Nebulae 

The condensation curves of a number of species are given for 
a range of temperatures and pressures which spans those of 
interTst in both the solar nebula (solid curve) and the jovian 
subnebula (dashed curve). The symbols on the solar nebula curve 
denote the conditions at the orbits of Mercury, Venus, Earth, 
Mars, and Jupiter at the time of maximum temperatures in the 
nebula. The symbols on the jovian subnebula adiabat denote 
similarly the condensation conditions at the orbits of Amalthea, 
lo, Europa, Ganymede, and Callisto. The jovian nebula is likely 
to be optically thin (and hence nearly isothermal) outside the 
orbit of Ganymede, as indicated by the horizontal dashed line. 
Condensation of ammonia hydrate in the jovian system is regarded • 
as marginal at best, and may well not be possible. 
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Of »rp.„tl„e for»«.o„ i, 3.105 t*o.. as large a. that l„ 

the nebula at the corresponding point in the condensation 
sequence. ,n addition, serpentiniaation in the subnebula become, 
theroodvnaoicallv possible at about boOK, coopered to 3 bOK in the 

lar nebula. it is therefore entirely plausible that serpentine 

may form readily by oas-solirt • 

y gas solid reactions in the subnebula, but that 

the Nineties ot its toro.tion in the solar nebula ere 

prohibitively Slow. The model strongly suggests that metallic 

iron Will be rare or absent in lo, and hence that the osygen 

fugacity will be buffered at k,- 

at far higner values than the Fe-FeO 

buffer. 

A recent petrological model of lo (Lewis lom » w 

tuewis, 1981 ) discusses the 

condition, necessary Tor sulTur magma generation and sulTur 
dlclde-driven volcanism. 1„ that discussion it is shown that the 
..quidsulTur magma and Its driving gas must be derived Trom. 
aource material devoid oT metallic iron and rich in o.idited 
It 1 , further maintained that ,ul f ur magmati so could not 

-ae on .0 until after loss of any primordial endowment of water 

carbon dioxide ^ 

This massive loss of volatiles is 

energetically plausible for lo hor-= t- 

tor lo because of its strong tidal 

heating, but is e>!tremely difficult for a cold 

^ massive planet 

xke Mars. It is for this reason that Mars and lo. even if they 

had Identical initial compositions, would have followed divergent 
evolutionary paths. 


The composition calculated tor Europe U, 


^5 slightly above the 
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water ice condensation curve, as required by our original 
constraints on the model. Thermal evolution of Europa as a closed 
system would cause decomposition of the phyl 1 osi 1 icate component, 
resulting in release of large quantities of water to the surface. 
This is a plausible source of the surface ice layer on Europa. 

Critique of the Model 

The first and most general test to which we may subject these 
results is to determine whether the absolute and relative masses 
of the Galilean satellites are acceptably in accord with 
observation. We can readily check the relative masses by 
augmenting the mass of each satellite up to full solar composition, 
by replenishing the missing volatiles in solar proportions. The 
material of each satellite is then spread over an annulus centered 
about that satellite's orbit and extending halfway to the orbit o'f 
each of its two neighboring satellites. The surface density of 
this material projected upon the symmetry plane of the nebula can 
then be plotted and compared to the solar nebula model. Figure 3 
presents the data on this mass distribution calculated as shown in 
Table 3. The slopes of several different power law relationships 
between the surface density and distance are shown for comparison. 
It can be seen that the actual mass distribution in the jovian 
system, when the actual densities and approximate compositions of 
the satellites are properly taken into account, fits a surface 
density which varies as the -2.3 power of jovicentric distance,, 
not the -1.6 power as derived from analogy with the solar system. 
Note that Lunine and Stevenson calculate a theoretical slope 
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Caption to Fig. 3 

Radial Distribution of Mass in Jupiter's System 
The nebular density in the protojovian subnebula, calculated 
in the manner described in the te>:t, is compared with four 
different functional relationships between the surface density 
and jovicentric distance. The mass distribution inferred in this 
manner is significantly steeper than that which was calculated 
theoretically by Lunine and Stevenson, or that found empirically 
for the solar nebula (see Table 4). 




log r (cm) 
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Table 3 

Minimum Mass of the Jovian Nebulat 


Satel 1 i te 

Mass 

Factor 

M( solar) 

r (annulus) 

A (annulus) 

Surface 

Density 

Amal thea 

<1 

*^240 

<240 

1.81 

2 / 

<1,000,000 

lo 

Q9. 1 

210 

18700 

4.22 

67 

2,790,000 

Europa 

48. 7 

. . 200 

9740 

6.71 

143 

680,000 

Ganymede 

149.0 

80 

11900 

10.71 

440 

270,000 

Cal 1 i sto 

106.5 

80 

8520 

18.83 

1285 

66,000 

J6,7, 10 each <1 

80 

<200 

116. 

70000 

<30 


24 10 20 ? 

♦ Masses in 10 g, radii in 10 cm, area in 10 cm 
3 -2 

density in 10 g.cm 


9 


surface 



r 
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of about -1.5 from their model of the dynamics of the nebula. We 
may attribute the discrepancy either ignorance of the fluid mass 
distribution in the jovian nebula or, more plausibly, to the very 
severe modification of the original distribution of solids by the 
effects of frictional drag between small particles and the very 
dense gas. A naive assessment of the effects of drag would surely 
attribute the most severe drag effects to condensates in the 
densest (innermost) portions of the nebula, and hence it would be 
reasonable to expect that friction would broaden and flatten the 
initial distribution of solids. This is the opposite to the 
effect actually observed. An alternative possibility is that the- 
rate of accretion is so phenomenally high in the inner portion of 
the subnebula that bodies rapidly accrete to sires large enough to 
be insensitive to drag effects. 

In order to make clear the basis for these comparisons between 
the solar nebula and the jovian subnebula, it is necessary to 
include a specific and explicit model for the distribution of 
matter in the solar nebula. Such a model, based on recent ideas 
regarding the compositions of the jovian planets, is given in 
Table 4. The mass factors needed to convert an observed solid 
body of approximately known composition into its solar-composition ' 
counterpart are shown in fable 5. 
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Table 4 

Minimum Maas of the Solar Nebula* 


Planet 

Mass 

Factor 

M(solar) 

r (annulus) 

A(annulus) 

Surface 
Densi ty 

Mercury 

3.3 

448 

1480 

0, 

1.82 

013. 

Venus 

48.7 

238 

11590 

0.83-1.29 

3.06 

3788. 

Earth 

59.8 

223 

13325 

1.29-1.89 

6.00 

2220.. 

Mars 

6.4 

216 

1382 

1.89-3.20 

20.95 

66. . 

<ast. ) 

(<0. 1) 

200 

<<20) 

3. 2-6.0 

80.9 

<0.3 

Jupiter 

14040. 

5 

70200 

6.0-11.0 

267. 

262. 

Saturn 

5695. 

12 

68380 

11.0-21.5 

1072. 

64. 

Uranus 

870. . 

44 

38230 

21.5-36.8 

2802. 

13.’7- 

Neptune 

1032. 

50 

51600 

36. 8—52. 0 

4240- 


Pluto 

0.01 

80 

<1 

52-70 

6900. 

4.6 

* Masses in 10 
—2 

26 13 

g, radii in 10 cm, area 

26 

in 10 cm , 

surface 


density in g,cm 
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Table 5 


Mass 

Factors -for 

Condensates 

in a Solar-Composition 

System 

Category 

No. of Moles Mass 

<Si=l) 

Mass Fraction M(i> Factor <=)/«) 
of Solar Mix 

Dry rock 


196.0 

0.0047248 

212. 

FeO 

0.6O0 

43.2 

0.0010414 


FeS 

0. <»00 

35.2 

0.0008485 


MgO 

1 . 060 

42.4 

0.0010221 


Si 02 

1 . 000 

60.0 

0.0014464 


CaO 

0.072 

4.2 

0.0001012 


A1203 

0.042 

4.3 

0.0001037 


Ni 

0.048 

2.6 

0.0000627 


Na20 

0.030 

1.9 

0.0000458 


All others — 

2.2 

0.0000530 


Ices 


577.0 

0.013909 

72. 

H2Q 

17-80 

320.4 

0.007724 


NH3 

3.74 

63.6 

0.001533 


CH4 

11.80 

188.6 

0.004546 


Ar 

0. 12 

4.4 

0.000106 


Gases 


40710. 

0.981365 

1.019 

H2 

15900.0 

31800. 

0.766580 


He 

2210.0 

8840. 

0.213099 


Ne 

3.4 

70. 

0.001687 

• 

Total 


41483. 

1.000000 

1.000 


A-i . 
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It is not presently possible to subject the compositional 
predictions of this model to very demanding tests, due to the fact 
that our present knowledge of the compositions of the Galilean 
satellites is limited to bulk densities and some (literally) 
superficial spectroscopic data. Several of the few available 
compositional constraints have already been mentioned in the 
previous pages, and no conflict between these limited data and the 
predictions of the present model is evident. 


Directions for Future Research 

The behavior of the- solids present in the dense jovian 
subnebula is of great interest. I have indicated above the 
existence of evidence for important radial transport of solids in 
the subnebula, presumably as a result of frictional dissipation 
caused by the velocity difference between solids and the embedding 
gas. This inward evolution of the solids, combined with the very 
short synodic periods of bodies in nearby orbits and the very high 
densities of dust and gas in the inner regions of the subnebula, 
suggest strongly that a detailed study of the dynamics of orbital 
evolution and accretion of solids in this system would be of great 
value. Interestingly, recent advances in the treatment of the 
motion and accretion of solids in the solar nebula both in the 
United States and Japan have provided a basis for such modeling. 

rther, the interesting guestion of the vapor and condensed phase 
transport mechanisms in the solar nebula has become the subject of 
joint effort by myself and R. G. Prinn (funded through grants to' 


a 
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M.I.T. ), and will surely be great importance in the jovian 

subnebula. In this study, the outward eddy transport of vapcrs 
across condensation fronts, partly balancea by the inward 
evolution of solid bodies due to viscous drag, leads to complex 
and markedlv nonuniform di str i but i cns of condensables in the 
system. Application of this approach to the jovian subnebula, 
once the principle has been worked out for the solar nebula, would** 


not be difficult. 
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Abstract 

oxkto" I"'l"<ll"8 AljCtj a, well „ heli<l„, 

SMcheiele I ^ ' l,'l ° “’"«<* «"l> theimodyiumlc end 

!™nttl ,r r 1 °^ o' 

<*se «l l ' e ' "" *" ”'" 'o >» loo ure to produce the 

.rved cloudt. Antlnron, me, ccoceleehly be tulflclentl, volatile, bu, the expected molecular 

o^ T.Ts“ tZ t »' '—“-v 

^ and Sb. Thus the arguments for a volatile-deficient origin for Venus based on the 

arseZn T" *” 

arsenic or antimony. ^ 



Soviet spacecraft have analysed cloud particles in the main Venus cloud layer by means of x- 
ray fluorescence (XRF) spectroscopy (1). These analyses claim detection of Cl as a cloud 
constituent, with S not detected. These jn^ results are in conflict with the weight of 
evidence from Earth-based studies (2) and from other spacecraft experiments (3), which 
strongiy suggest that the dominant cloud constituent on Venus is H2SO4 droplets. Satisfactory 
photochemical models for the production of an H2SO4 aerosol from geochemically plausible 
primary gases, including COS, H2S and SO2 are available (4). Unfortunately, direct evidence 
regarding the abundances of these species In the lower troposphere is lacking: copious 

production of COS, etc. by reactions between sulfuric acid and the inlet system of the mass 
spectrometer on the Pioneer Venus large probe leads to masking of the atmospheric sulfur gases, 
and abundances up to 100 ppm are possible for COS and H2S (5). 

The source of a chlorine-bearing aerosol is less obvious. We have previously pointed out the 
high volatility of halides and sulfides of mercury, arsenic and antimony (6). and have shown that 
the terrestrial crustal abundance of even the rarest of these elements. Hg. would suffice to 
produce substantial masses of halide cloud condensates on Venus: mercury is so volatile at the 
surface temperature of Venus that It would reside almost completely in the atmosphere. 


More recently, in the context of a model for the formation of the planets in the presence of 
a steep gradient outward from the proto-Sun, we have favored compositional models in which 
the volatile elements are severely depleted in the accreting Venus relative to Earth (7). The 
failure of the 1978 Pioneer Venus mission to detect even a trace of Hg in the lower atmosphere 
(5) strongly implies that Venus Is deficient in Hg relative to Earth. The severe depletion of the 
most abundant terrestrial volatile, HjO, on Venus is well known, and is variously attributed to 
either a lack of water in pre-planetary solids at the orbit of Venus (7,8) or to massive loss of 
oceans worth of H2 after differentiation and outgassing of the planet (9). Such a late loss 
mechanism could not, of course, deplete mercury while leaving vast amounts of the lighter and 
more volatile species N2 and COj. Thus the observed severe depletion of Hg on Venus is more 


convincing evidence for a volatile-poor high-temperature origin of Venus than is the depletion 
of water. Other moderately volatile elements, such as arsenic, antimony, bismuth and 
germanium are also potential indicators of the overall volatile content of Venus. In addition, 
such species, if present in the hot lower atmosphere would condense at intermediate altitudes to 
form solid halides, sulfides, and oxide cloud particles. 


In the past year, Krasnopolsky and Parshev (10) have suggested AI2CI6 as the major cloud 
layer constituent, and Barsukov et al (11) have suggested arsenic and antimony halides and 
oxides. Mole fractions of 0.1 to 1 ppm of condensible gases are required to provide the 
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obseK^ed cloud density (12). Are these species plausible cloud constituents? Can useful limits 

on the abundances of the volatile elements As and Sb be derived from atmospheric 
measurements? 

It is a simple matter to estimate the volatility of aluminum compounds at the mean surface 
conditions of Venus. Consider co-existing anorthite (CaAl2Si208), quartz (Si02) and calcite 
(CaC03) 730 K and 95 bars in contact with the atmosphere: 

CaAl2Si208 ♦ C02(g) ♦ 6HCI(g) = CaC03 ♦ 2Si02 * 3H20(g) ♦ 2AICI3 (1) 

The equilibrium constant for this reaction at this temperature is 10-26.1. Assuming the 

spectroscopic HCI abundance in the clouds, about 1 ppm (13), and the most typical water vapor 

abundance figures for the lower atmosphere, about 100 ppm (14), an AICI3 partial pressure of 

10 21.0 b is calculated. With an extreme effort to bias the equilibrium In favor of AICI3 

production (raising HCI to 10 ppm and lowering H2O to 10 ppm), we can force the AlCI, pressure 

only as^ high as 10 16.5 b. This is still 10 orders of magnitude too small to provide detectable 

amounts of condensate.- Note that the presence of granitic rocks on Venus, with anorthite and ' 

quartz as ccxnmon primary minerals, is expected both from considerati-Dn of the atmospheric 

composition (15) -nd from surface passive gamma-ray spectroscopy (16). The surface of Venus 

lies, as accurately as can be determined, precisely on the calcite-quartz-wollastonite (CaSiO,) - 
buffer: '3' 


CaC03 * 5*02 = CaSi03 * 002(g) (2) 

and calcite is therefore also a plausible surface mineral. Since weathering reactions at higher 
altitudes on Venus will preferentially tend to mobilize a fine Ca-rich dust, which can be 

transported readily by winds to the hot lowlands, this buffer may in fact not be difficult to 
establish on Venus (17). 

Thermodynamic treatment of arsenic and antimony volatilization can be carried out without 
a precise a priori_ knowledge of the minerals formed by these elements on the surface (18,19). 
We shall calculate the partial pressures of a number of As and Sb gases at the surface as a 
function of the activities of As and Sb: an activity of one means the pure element is present on 
the surface; and activity of 10'^ means that the pressure of the monatomic vapor of that 
element is 10-» times its abundance at saturation. We can then assess the stability of possible 

surface minerals containing these elements. Hg. 1 presents the results for arsenic (20). For As 

activities greater than.about 10'2 the dominant gas is AS4 with AS4O8 second. Based on the 



Captions of Figures 


Fig. 1 . Partial pressures of arsenic gases as a function of elemental arsenic activity at 
the Venus surface. Limitations of the arsenic activity imposed by precipitation of liquids 
of As 2$3 (orpiment) and AS 4 S 4 (realgar) composition are indicated by the vertical dashed 
lines. Heavy dots indicate the upper limits on the AS 4 , AS 4 O 6 , and AsS partial pressures. 
A condensible species with a partial pressure near 10*4 bars is needed to provide the 
observed cloud density. The maximum total mole fraction of arsenic gases is i 0.1 ppm. 

Fig. 2 . Partial pressure of antimony gases as a function of elemental antimony activity. 
The Sb 40 g(valentinite) and Sb 2 S 3 (stibnite) precipitation points are Indicated by the 
vertical dashed lines. Firm upper limits on the antimony gas pressures are Indicated by 
the dots. The maximum total mole fraction of antimony gases is about 0.03 ppm. An SbS 
gas abundance of 0.3 ppm is needed to make clouds. 
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recfot thermochemical data of Johnson et al (18), we calculate that liquid AS2S3 will precipitate 
If t‘.e elemental arsenic activity 1$ greater than 0.13, thus ruling out higher As activities. This 
In turn places firm upper limits on the pressures of As-bearing gases as indicated in Fij;. - 

The results of a similar calculation for antimony are given in Fig. 2. Pre ip.:ition of 
Sb406(s) occurs for an elemental antimony activity of 0.20. Thus for all possible Sb activities, 
SbS is the dominant gas. The best available data on Sb2S3(stibnite) from Johnson's work (18) 
Indicate that Sb2S3{s) precipitates at an antimony activity of 10*2*3. Therefore we expect that 
stibnite precipitation on the Venus surface will regulate the Sb gas phase abundance. 

We conclude that the mole fraction of all arsenic gases is below 10*2 probably making these 
species too rare to account for the clouds no matter what species condenses. We expect that 
the total mole fraction of antimony gases will be lower than for arsenic, but that the most 
Stable gas is SbS, not a halide. Partial pressures as high as 10*^ bars are conceivable for SbS 
and cannot be confidently ruled out. 

Observational constraints on the abundance of arsenic in the Venus atmosphere even down to 
the level of 0.1 parts per million are therefore not sufficient to test whether arsenic is, like 
water and mercury, depleted in Venus relative to Earth: the stability of arsenic sulfides is great 
enough to preclude a larger abundance of gaseous arsenic compounds irrespective of the crustal 
abundance of As. The same may be true of antimony, since the mineral Sb2S3 seems to have 
low_ enough volatility to hide Sb in the lithosphere. We have also briefly considered bismuth, and 
find that Bi2S3 is so stable that the most abundant Bi bearing gas, BiS, should have a mole 

fraction below 

In any event, we call into serious question the geochemical plausibility of al^the species so 
far suggested as sources of chlorine-bearing clouds. We suggest that either the chlorine 
compound is a species which has not been considered, or the XRF data used to deduce the 
presence of chlorine may be in error (21). 
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Abstract 

A set o-f e-f-ficient programs for calculation of condensation 
behavior in a system with either solar or carbon-rich elemental 
composition has been developed to treat the course of condensation 
at very low pressures. Condensation processes of astrophysical 
interest, such as those occurring in the atmospheres of C stars 
and in nova and supernova envelopes, can be treated by these 
programs. These pro, ’.rams have been applied to the problem of 
condensation at very low pressures, where intergrain .ollisions 
are rare and nucleation sites for silicates are found only on the 
surfaces of previously condensed grains. This "isolated grain" 
case gives results which d .-erge in interesting ways from both the. 
equilibrium condensation ano the nonhomogeneous accretion models. 

Extension of the progr. ms to a set of 20 elements is under way, 
and a full set of publishable results for the do '.densati on process 
in the above cases of astrophysical importance will soon be 
available. This project is also producing calculation procedures 
which make it possible to carry out complex thermochemical 
equilibrium calculations using a small computer. 

The minerals produced in the stellar and nova-related processes 
under study, including carriers of important volatile elements 
such as carbon and. nitrogen, are candidates for accretion into 
meteorite parent bodies and planets, and may still be discernible 
in the enstatite chondrites. 
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Introduction 

Considerable work has been done during the past year on the 
condensation behavior of a solar-composition gas in cool stellar 
environments and in carbon-rich astrophysical environments such as 
C nova shells. These environments are the original sites of 
condensation of the solids in interstellar clouds and preplanetary 
nebulae, and hence are the ultimate sources of the solids present 
in the solar nebula. The first phase of this project concentrated 
on developing FORTRAN programs to treat the gas and condensation 
chemistry of solar material, with the novel constraint that 
reactions between two or more dissimilar grains and the cooling 
gas were omitted. This corresponds to the physical condition 
that grains formed in astrophysical settings generally do so at 
extremely low pressures, usually less than 0.01 dynes per square 
centimeter. Thus nucleation sites are at a premium, grains almost 
never collide, and condensation takes place at unusually low 
temperatures. Therefore fresh condensates tend to coat old 
grains, and diffusion is very slow. Disequilibrium between the 
solids in the grain interiors and the gas is enhanced bv the high 
diffusion barrier, but also by the short dynamic expansion and 
cooling times inherent in nova and supernova ejecta and in carbon 
star "smoke". Some of the early condensates, especially in 
system with an excess of carbon over oxygen, are excellent 

Candidas, s as carriers of volatile elements, especially carbon 
and nitroqen. 

The puczlinq departures from equilibrium seen in certain 
classes of meteorites, especially the carbonaceous chondrites and 
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the enstatite chondrites, point strongly to non-nebular processes. 
In the case of the carbonaceous chondrites, secondary alteration 
processes in a parent body have become widely accepted. Secondary 
alteration processes have also been proposed to explain the 
strange coexistence of very highly reduced high-temperature 
silicates with abundant volatiles and low-temperature sulfides in 
the enstatite chondrites, but this scenario has not been generally 
accepted. Both classes, in fact, contain carbon and other highly 
reduced phases which suggest origins in astrophysical settings 
different from the classical simple-minded picture of a solar 
composition nebula. As mentioned in Section I, we have already 
begun to look at gas and dust transport processes in a dynamical 
model of the solar nebula to see whether such transport may 
contribute to the stability of the odd phases in chondrites. The 
present project deals with the possibility that these odd grains 
may have been inherited from a presolar interstellar cloud without 
extensive chemical alteration. 

Program Development and Computation Strategy 

During the initial phase of this project, calculations were 
carried out from a terminal leased by the company, using the 
Multics system at the MIT computer center. An arrangement was 
made by which c. p. u. time would not be charged on this project. 
That arrangement was withdrawn by MIT last summer, so that 
after-hours work was necessary to realise reasonable rates. With 
the availability of the new IBM Personal Computer, it- became 
possible to remove all computation from large, expensive systems 
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and thereby do away with connect-time and cpu-time charges. This 
was so attractive a possibility that we immediately translated 
some of the fundamental computation routines into Basic for use on 
the IBM. This conversion was carried out very efficiently, and we 
are now running in Basic. A FORTRAN compiler has been announced 
for the IBM, but is not yet commercially available. 

The present strategy employs gas equilibration routines which 
can be stepped in both temperature and pressure, and which 
calculate the condensation threshholds of a number of solid 
phases. The program is run interactively, so that the appearance 
of a new condensate can be screened for its stability under the 
rules of the procedure < isolated-grain behavior), and accordingly 
accepted or rejected. 

The computer programs used ir, this project differ substantially 
from the f ree-energy-minimization routines ' which have become 
widely used. They have the advantage that they are readily 
adaptable to use on small computers without using excessive core, 
and that un;'er almost ary circumstances they compile and execute 
in a fraction of the time. The procedure is an iterative 
improvement .:.f elemental activities (initially estimated by a kind 
of programme chemical intuition) in which the magnitude of the 
corrections io the activities is determined by the comparison of 
the total computed abundances of the compounds of each element to 
the initially specified elemental abundance of that element. A 
pressure-temperature point can be solved in a system of 16 
elements within one to tw-> minutes by the IBM Personal Computer, 
and this performance, can be realized with programs which have not 


been purged oi several obviously inefficient (but reliable ) 
procedures. The Basic program is now being expanded to 20 
elements, which is the size of the largest FORTRAN program of this 
general type so far developed (the program T0P20, developed by S. 
S. Barshay while he was a student of mine at MIT, which could only 
be run on CDC6600 or larger computers!). The problem of exponent 
overflow and underflow, a perennial issue with IBM computers, can 
be completely and efficiently bypassed by a new scaling procedure 
which ws h^vG dsvelopGci i^his 

We have already reported briefly on some of the features of the 
isolated-gram condensation process in our renewal proposal, and 
we include here a summary of the more interesting features of this 

process which set it apart from the more familiar equilibrium 
condensation process. 

Results to Date 

The general features of the isolated-grain condensation process 
which we have documented to date include important differences in 
the behavior of iron and the alkali metals. The nature of the 
condensation sequence can best be appreciated by comparing the 
results to the more conventional equilibrium condensation and 
nonhcmoqeneous accretion schemes. 

The first condensates in all three schemes are refractory oxide 
minerals rich in aluminum, calcium, and titanium and poor in 
silicon. These include corundum, spinel, perovskite, anorthite 
and melilite sol id solution. The next condensate in all three 
schemes is metallic iron-nickel alloy. The third condensate is 



magnesium silicate. At nigh total pressures, forsterite is the 
first magnesium silicate to appear, followed by alteration to 
enstatite by reaction with si 1 icon-bearing vapors at slightly 
lower temperatures- At pressures of 0.1 dyne per square cm, the 
condensation temperature of iron has dropped below that of 
enstatite, so that enstatite slightly precedes metal as a 
condensate. The equilibrium condensation model predicts that all 
of the silicon will be used up in forming silicates, with none 
left over to form quartz. However, both the nonhomogeneous 
accretion model and the i sol ated-grain model effectively paint 
lower temperature condensates on high temperature condensates. 
Thus the former process rapidly accretes forsterite, then 
enstatite, then quartz onto growing pi anetesi mal s, while the 
latter similarly zones individual grains. The alkali feldspars 
require intimate contact between gas, anorthite, and a silica-rich 
phase such as quartz. This is prevented in the nonhomogeneous 
accretion scenario by early accretion and burial of anorthite, and 
in the isolated-grain model by the coating of anorthite grains 
(the main source of aluminum) by impervious layers of magnesium' 
silicates and silica. 

Note that, up to this point, the likely mineral products of the 
nonhomogeneous accretion and isolated grain models are essentially 
identical. The difference lies in the way that these mineral 

components would be distributed in a p 1 anetesi mal : the 

nonhomogeneousl y accreted body will of course be strongly layered, 
with a core rich in refractories and metal, and a silica-rich 
surface?. A body accreted froin the products of i sol ated-qr ai n 
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condensation would be chemically homogeneous on a distance scale 
of millimeters. An important point is that the equilibrium 
condensation products are quite distinct -from those in either of 
these two schemes, and more closely resemble ordinary chondrites, 
in that they contain two feldspars, no free silica, and have 
significant traces of FeO in the olivine and pyroxene. 

At this point the nonhomogeneous accretion and isolated grain 
condensation sequences diverge very significantly. Metallic iron 
IS of course buried and chemically isolated from the gas in the 
former model, but present as metallic grains in the latter. Thus, 
at the appearance temperature of troilite (6aOK) , it is possible, 
for the isolated-grain process to make FeS unless silicates have 
nucleated on the metal grains. In view of the presence of 
silicate grains at all temperatures at which ' metal is condensed, 
it 13 plausible that nucleation of 1 ower-temperature silicates 

will strongly favor these silicates, and thus not coat metal 
grains. 


Phosphorous, like sulfur, should react in the isolated grain 
case with metal grains, but, due to the limited ability of these 
grains to absorb phosphorous and retain it at lower temperatures, 
phophates should be ultimate products in all three schemes IF 
calcium were chemically available to form phosphates. 

The issue of the behavior of phosphorous, along with the 


co/riple>: 


xses involving the alkali sulfides 


and halides, is still under stud'/ 
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The following program listing is illustrative of the present 
level of development of the Basic translations and adaptations of 
our original FORTRAN programs, which were suitable only for use on 
very large computers. 


10 

11 

20 

26 

27 

28 
29 

40 

41 

42 
50 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 


READ BH, BHE, BQ. BC, BN, BSU, BNE, BAR 
READ BSI, PFE, BMG, BAL, BCA 
FOR NT « 1 TO 10 

READ T, KH, KH20, KCO, KCH4, KCQ2 
READ T, KOH, KO, KC, KH2S, KHS, KN 
READ T, KCN, KHCN, KNH3, KS02, ‘ KS 
READ T, KSIO, KSI, KSIH4, KFE, KMQ 
READ T, KCA, KAL, KALO 

READ T, KSI02, KMGSI03, KMGO, KMG2SI04 
READ T, KFES, KFESI03, KFE2SI04, KFEO 
READ T, KAL203, KCAO, KSPINEL, KAN. 25 
FOR NP = 1 TO 10 
F = 10! ''(NP-10! ) 

ABH = F*BH 
ABHE = F*BHE 
ABD = F*BQ 

ABC =» F*BC . \ . 

ABN = F*BN 
ABSU = F*BSU 
ABNE = F*BNE 
ABSI = F*BSI 
ABFE = F*BFE 
ABMG = F*BMG 
ABAL = F*BAL 
ABCA = F*BCA 
ABAR = F*BAR 


8o REM INITIAL GUESS OF ELEMENTAL ACTIVITIES 
31 RH = -(KH/4) +SQR <KH''2+Q*ABH) 

82 RN = -(KN/4)+SDR(KN-'-2+S*ABN> 

83 RO = (2*AB0) / (3*KH20*RH-'2) 

84 PH2 = RH'-2 

05 GR = ABC/ ( (KC0*R0) (KCH4*PH2"2) ) 

86 RS = AE<SU/ < (KSQ2*R0''2> + (KH23*PH2) ) 

90 ASI = ABSI / <KSI0*R0) 

92 AMG = ABMG*KMG 

93 AAL = ABAL/KAL 

94 ACA =• ABCA /KCA 
99 INDEX = O 


200 FOR ITER = 1 TO 150 

201 INDEX = INDEX-t-1 

202 PH = KHKRH 

203 PH2 = RH -2 

204 PHE = ABHE 



211 PCH4 = KCH4*GR*PH2-^2 
212.PC02 = KCQ2*GR*RQ''2 
.213 PC = KC*GR 

220 POH = KOH*RH*RO 

221 PH20 = KH20*PH2*R0 

222 PO = ka*RQ 

223 P02 » RQ''2 

230 PCN » KCN*GR*RN 

231 .PHCN =* KHCN*GR*RH*RN 

232 PNH3 a KNH3!KRN*PH2*RH 

233 PN2 = RN'-'2 

234 PN = KN*RN 

240 PH2S = KH2Si«RS*PH2 

241 PSQ2 = KS02*RS*P02 

242 PHS = KHS*RH*RS 

243 PS = KS*RS 

244 PS2 = RS''2 

250 PSIO = KSIO*ASI*RO 

251 PSl = KSI*ASI 

252 PSIH4 = KSIH4*ASI*PH2-''2 
250 AFE = ABFE/KFE 

259 IF AFE>1! THEN 260 ELSE 264 

260 PFE => KFE 
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261 PMETFE = ABFE - PFE 

262 AFE = 1 ! 

263 GOTO 270 

264 PFE = ABFE 

265 PMETFE = 0 ! 

270 PMG = ABMG 
280 PCA = ACA^KCA 

201 PAL = AAL*KAL 

202 PALO = KALO*AAL*RO 

300 REM CALCULATE ACTIVITIES OF SOLIDS 

301 ASI02 = KSID2*ASI*P02 

302 AMGSI03 = KMGSI03*AMG*ASI *P02*R0 

303 AMGO = KMG0*AMG*R0 . 

304 AMG2SIQ4 = KMG2SI04*AMG'''2!KASI*P02'^2 

305 AFESI03 = KFESI03*AFE*ASI *P02*R0 

306 AFE2SI04 = KFE2SI04)((AFE-'2*ASI *P02 '2 

307 AFES = AFE*KFES*RS 
303 AFEO = KFE0*AFE*R0 

320 ACAO = KCAO*ACA^;RO 

321 AAL203 = KAL203*AAL''2*P02*R0 

322 ROOTAN = KAN. 25*P02*SQR ( AAL^ASI *SQR ( ACA) ) 

323 AAN = R00TAN-'4 

324 ASPINEL = KSPIMEL*AMG*AAL'-2*P02''2 
40C REM CALCULATE ELEMENTAL SUMS 

401 SH = PH+2*PH2+2*PH20+4*PCH4+3*PMH3+P0H+PHS+PHCN+2*PH2S 

402 SC = PC0+PC02+PCH4+PC+PHCN+PCN 

403 SO = PH2Q+2*P02-*-2)KPC02+PC0+P0H+2*PS02«-PSIQ+P0+PAL0 

404 SN = PN+2*PN2+PNH3+PCN+PHCN 


405 SS = PS+PHS>PS02+PH2S+2*PS2 


410 SSI = PSI0+PSI+PSIH4 

411 bMG = PMG 

412 SFE = PFE + PMETFE 

413 SAL - PAL+PALD 


414 SCA = PCA 

4 1 9 p = PHI’ < PH+PH20+PC0 ♦ PC02-i PCH4+PM2+PNH3+PaH+P = lO+PFE *-PMG+PH2S+PHE+P0+PNE+PHS+ 

PAL+PCA+PALO+PAR . . 

420 EH = ABS <SH-ABH) /ABH • •• 

421 EQ = ABS <SO-ABO) /ABO 

422 EC = ABS (SC-ABC) /ABC 

423 EN = ABS 'iSN-ABN) /ABM 

<’ c:q— AC vGi n /amCM 



431 EFE =* rtBS(SF£-ttEiE£)/ABFE 

432 EMG = ADS(SMG~ABMG)/ABMG 

433 EAL = ABS(SAL-ABAL>/ABAL 46 

434 EGA => ABS<SCA-ABCA) /ABCA 

440 IF EH<.001 THEN IF E0<.001 THEN IF EC<.001 THEN IF EN<.001 THEN IF ES<.001 T 
HEN IF ESK.OOl THEN IF EFE<.001 THEN IF EMG<-001 THEN IF EAL<.001 THEN IF ECA< • 
001 THEN GOTO 501 ELSE GOTO 449 

449 REM COMPARE ELEMENTAL SUMS TO ABUNDANCES 

450 IF INDEX - 1 THEN 451 ELSE 455 

451 AFE = AFE*(ABFE/SFE) 

452 RH => RH*SQR(ABH/SH) 

453 AAL => AAL* (ABAL/SAL) ORIGiNAL PAGE !S 

454 GOTO 500 OF POOR QUALITY! 

455 IF INDEX = 2 THEN 456 ELSE 460 

456 RN =» RN*SQR<ABN/SN) 

457 RO RO* (ABO/SO) 

459 GOTO 500 

460 IF INDEX = 3 THEN 461 ELSE 4to5 

461 AMG = AMG*(ABMG/SMG) 

462 GR == GR* (ABC/SC) 

463 ACA = ACA*(ABCA/SCA) 

464 GOTO 500 

465 RS = RS*<ABSU/SS) 

466 ASI - ASI* (AB3I/SSI) 

496 INDEX = 0! 

500 NEXT ITER 

501 REM OUTPUT RESULTS OF ITERATION 

502 BEEP 

505 PRINT "TEMPERATURE =* ";T, "PRESSURE - ";P, ITER; " ITERATIONS" 

507 PRINT " PRESSURES OF GASES IN BARS" 


508 

PRINT 

"H =";PH; " 

H2. 

=";PH2; " 

HE =";PHE;" 

NE =";PNE;" AR =»";PAR 

509 

PRINT 

"0 =";P0; " 

02 

=";P02; " 

OH =";POH;" 

H20 =";PH20 

510 

PRINT 

"N =";PN?" 

N2 

=";PN2; " 

NH3 =";PNH3; 

" CN =";PCN 

511 

PRINT 

"C =";PC;" 

CO 

=";PCO; " 

C02 =";PCQ2; 

1. =";PCH4;"HCN =";PHCN;" S02 


;PS02;" ALO =";PALO 

512 PRINT "S =";PS;" S2.=";PS2;" HS =";PHS;" H2S ; PH2S 

513 PRINT "FE =";PFE;" MG =";PMG;" CA =";PCA;" AL =";PAL 

514 PRINT "SIO =";P3I0;"- SI ==";PSI;" SIH4 =";PSIH4 

516 PRINT " ACTIVITIES OF CONDENSATES" 

517 PRINT "QZ =";ASI02;" MGO =";AMGO;" EN =";AMGSI03;" FO ; AMG2SI04; " * FE 
;AFE;" CAO =";ACAO;" COR =";AAL203 

5ia PRINT "FED =";AFEO;" FS ; AFESIQ3; " FA ; AFE2SIQ4; " FES =";AFES;"SP =";ASP 
INEL;" AN =";AAN 

520 PRINT " INPUT ABUNDANCES OF ELEMENTS" 

521 PRINT "H =";ABH;" 0 =";ABQ;“ C =";ABC;" N =";ABN;" S =";ABSU 

522 PRINT "FE =";ABFE;" SI ="ABSI;" MG =";ABMG;" AL =»";ABAL;" CA =";ABCA 

523 PRINT " CALCULATED ABUNDANCES OF ELEMENTS" 

524 PRINT "H =";SH;" 0 =";S0;" C =";SC;" N =";GN;“ S «";5S;"CA =’;SCA 

525 PRINT "FE ==";5FE;" SI =";SSI;" MG =":3MG;" AL =";SAL 

545 INPUT GOON 

546 IF GOON = 1 THEN 590 ELSE 650 
590 NEXT NP 

600 NEXT NT * 

650 END 



Directions for Future Work 

Our efforts during the first year of this contract have been 
concentrated on writing and translating the computer programs to 
carry out the task, and only recently have we reached the point 
where publishable scientific calculations may be carried out. It 
is our expectation that the crucial scientific questions involving 
the details of the i sol ated—grain condensation process and the 
comparison of this process with previously explored scenarios can* 
all be addressed directly within a month or two. 

We therefore plan to examine the condensation processes for 
minerals in both solar — composition and carbon-rich systems at very 
low pressures, in accordance with our original proposal. It will 
thwn be a simple matter to compare these results to the observed 
mineralogy of meteorites to see whether such materials have* 
persisted to the present. We are especially interested in 
identifying mineral carriers for the volatile elements, especially 
carbon, nitrogen and the halogens because of their 
possible importance in contributing atmophile elements to the 
accreting terrestrial planets, including the possibility that some 
of the carriers of rare gases in meteorites may have originally 
condensed in astrophysi cal , rather than solar nebular, settings. 
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